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On  the  Non-existence  of  Limiting  Lines 
in  Transonic  Flows 

By  C.  S.  MORAWETZ  and  I.  I.  KOLOONER 

We  consider  flows  past  airfoils  which  are  irrotationai  and  subsonic  for  low 
Mach  numbers  at  large  distances  from  the  airfoil.  If  the  Mach  number  at 
infinity  is  gradually  increased,  at  some  particular  value  a small  supersonic 
region  appears  in  the  flow  next  to  the  airfoil.  The  flow  appears  to  remain  without 
shocks.  If  the  Mach  number  at  infinity  is  increased  still  more,  a definite  observ- 
able shock  wave  appears.  The  question  arises  as  to  why  this  shock  wave  appears. 

It  is  possible  and  even  quite  probable,  that  no  mixed  flows  without  shocks 
exist  as  transition  stages  and  that  in  general  as  soon  as  the  supersonic  region 
appears  there  is  a shock  wave  which  is  at  first  so  weak  as  to  be  unobserved. 

However,  if  these  continuous  mixed  flows  with  varying  Mach  number  do 
exist,  an  explanation  of  why  the  continuous  flow  breaks  down  is  needed,  and 
it  has  been  suggested  by  various  authors  that  a “limiting  line”  appears  in  the 
flow.  Tollmien,  Ringleb,  von  KArmAn,  and  Tsien  have  observed  that,  if  flows 
and  their  corresponding  profiles  are  constructed  using  solutions  in  the  hodo- 
graph  plane  which  depend  continuously  on  some  Mach  number,  there  is  always 
a critical  Mach  number  where  the  mapping  of  the  hodograph  plane  into  the 
physical  plane  breaks  down.  In  fact,  the  image  in  the  physical  plane  has  a fold 
in  the  supersonic  region  whose  edge  is  known  as  the  “limiting  line.”  The 
Jacobian  cf  the  transformation  from  the  hodograph  plane  to  the  physical  plane 
vanishes  along  the  limiting  line. 

Tsien  [1]  and  von  KArmAn  have  proposed  that  continuous  flow  past  a 
fixed  airfoil  also  breaks  down  because  of  the  appearance  of  a limiting  line. 
Nikolskii  and  Taganov  (2]  have  shown  that  such  a limiting  line  would  have  to 
8 tart  on  the  sonic  curve. 

It  was  finally  shown  by  Fiiedrichs  [3]  that  limiting  lines  cannot  appear 
anywhere  in  analytic  flows  which  depend  continuously  on  the  entrance  Mach 
number  and  that  therefore  the  breakdown  of  potential  flow  must  be  due  to 
other  causes.  (Friedrichs’  proof  was  challenged  in  a controversial  review  by 
Tsien  [4].)  Manwell  [5]  has  shortened  the  proof  considerably  and  eliminated 
the  difficult  lemma  supplied  by  Flanders  in  [3].  The  proof  presented  here 
dispenses  with  the  condition  of  analyticity  and  requires  only  continuous  second 
derivatives  of  the  stream  function.  In  addition  we  shall  show,  as  in  (5),  that 
in  the  construction  of  flows  past  continuously  changing  profiles,  [6,71  etc.,  the 
incidence  of  a limiting  line  corresponds  to  a profile  of  infinite  curvature. 

These  two  results  are  formulated  in  two  theorems. 
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Theorem  1.  A limiting  line  cannot  appear  in  a plane  continuous  flow  past 
an  airfoil  of  bounded  curvature  if  the  flow  depends  continuously  on  the  Mach  number 
at  infinity  and  has  a bounded  supersonic  region. 

Theorem  2.  If  a set  of  flows  which  depend  cordinuously  on  a parameter  is 
constructed  in  the  hodograph  plane  together  with  the  corresponding  set  of  profiles 
and  if,  for  some  critical  value  of  the  parameter,  a limiting  line  appears,  then  the 
corresponding  profile  has  infinite  curvature  at  some  point. 

We  make  the  following  assumptions:’ 

(a)  A solution  in  the  hodograph  plane  of  the  equations  of  an  isentropic, 
irrotational,  plane  flow  is  given ; i.e.,  we  have  a potential  tunction  <j>(q,6 ) and  a 
stream  function  f(q,6)  which  satisfy  the  equations 

<t>$  “ 

U) 

= P V‘e  *(?*  — c')i, 

in  a closed  simply  connected  region  D of  the  7,6-plane.  Here  p is  a known  func- 
tion of  q obtained  from  the  Bernoulli  theorem  for  steady  isentropic  irrotational 
flow,  and  c*  — dp/dp,  where  p — p(p)  is  the  equation  of  state  relating  pressure 
to  density. 

(h)  f,i  , 4>u  , exist  and  are  bounded  in  D. 

(c)  The  boundary  C of  D intersec  ts  the  sonic  line  q — c « cm  in  exactly 
two  points. 

(d)  p > 0 throughout  D. 

(e)  The  equation  of  state  p — p(p)  is  such  that  d*p/dp*  > 0. 

A solution  + of  equations  (1)  gives  rise  to  a flow  in  the  physical  plane 
provided  that  the  Jacobian  J of  the  transformation 


x - x(q,6) 


y = y(qfi) 

does  not  vanish.  J is  given  by 


(2) 


d(e,q)  d(q,6)  \d(r,y)/ 


“ -j-  - (7*/c*  — i)?'VD  * ?"*($•  - P'VV  - 

P 2 


A limiting  line  occurs  if  J vanishes  along  a curve.  This  cannot  happen  in 
the  subsonic  region,  see  [2],  [3],  and  [8].  Thus  a limiting  line  can  occur  only  in 
the  supersonic  region  or  on  the  sonic  curve. 

In  connection  with  J we  now  prove  the  fundamental 


'Standard  notation  for  the  flow  quantities  is  need. 
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Lemma : 

If  / > 0 on  C then  J > 0 in  D for  q > e,. 

Proof: 

We  introduce  the  new  dependent  variables 

U *»  9*  — tyt 

(3) 

V = <t>,  + k\f/t 

where 

(4) 

k * = p_Io“*(g*  — c*). 

With  this  change  of  variables,  J becomes,  by  (2) 

(6)  J « q~'UV\ 

hence  J «*  0 if  and  only  if  U — 0 or  V -»  0 since  by  assumption  (d)  both  U and 
V are  bounded  in  D. 

From  the  differential  equations  (1)  we  find  differential  equations  for  U and 
V.  First  we  introduce  the  characteristic  variables  a and  0 by  the  equations 


2 da  = d*  — dO 
9 

(6) 

2 d0  = &dq  + dO. 

q 

Then  we  have 


(7) 


± D jlA  _ A 

da  pk  dq  09 


0 ^ 0_  , 0_ 
08  " ok  da  + 00 

and  from  differentiating  (3)  we  find,  using  (1), 


Uf  =*  p k qU>t,  — — k,4/t)  4*  4>ti  — k^n 

* -P~lk~lqk,t,  , 

or  by  (3) 

(8)  U,  - hp-'k-'qk'iU  - V ). 

Similarly, 

(9)  Vm  = \p-'k-'qk,(V  - U). 

But,  by  (7),  since  k is  a function  of  q alone,  we  have 


fc,  = pq~lkkm  = pq'kkt  . 
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Substituting  in  (8)  and  (0),  we  thus  obtain  the  equations 


(10) 


U , - \k-%U  - -hp-'h-'qk.V 

V . - hk-'kaV  « —\o~'k~iqkrU 


which  may  be  rewritten  as 

(ID 

where,  by  (4),  (10), 


(fc-'U*), 


(*"'  V 


A 

/• 


-BVV 

-BUV 


Thus  B > 0 for  q > e„  > c since  p,  < 0 by  Bernoulli’s  theorem  and  dc2/dp  — 
d*p/dp*  > 0 by  assumption  (e). 

Let  us  suppose  first  that  J vanishes  or  is  negative  somewhere  in  the  super- 
sonic region  q > c„  > c of  D.  Since  J is  continuous  and  D is  closed,  there  is 
a maximum  value  of  q > c„  for  which  J — 0.  Hence  by  (6)  and  (4),  there  is  a 
maximum  value  r = a*  + d*  of  a 4-  d for  which  J *=  0,  and,  for  a + d > rt 


J > 0,  since  J > 0 on  C.  By  (5),  either  U — 0 or  V **  0 for  a ■=  a*,  d — d*. 
Suppose  U = 0 at  this  point.  Then  k~lU*  = 0 at  (a*, ft*)  but  k~'U 7 > 0 for 
a + d > t.  Hence,  for  a — a*,  d > d*>  A:-* f/*  > 0 and  therefore,  by  the 
mean  value  theorem,  (kTlU*)t  > 0 for  a = a*  and  some  value  d of  d where 
d > d*«  But  then  by  (11)  since  B>  0 we  have  UV  < 0 at  (a*,  d),  that  is,  by 
(5),  J < 0 at  (a*,d).  Since  a*  + d > r this  contradicts  / >0fora-fd>T 
and  hence  U cannot  vanish  for  q > c#.  Similarly  V cannot  vanish,  and  thus  J 
cannot  vanish  for  q > c*  > c. 
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Next  fippoec  tha*  j - 9 at  a point  P where  q{P)  — c#,  and  let  a(P), 
0 (P)  be  th  • coordinates  o.  this  point  in  the  characteristic  coordinate  system. 
By  (11),  sir  ne  7 > 0 an-;  hence  UV  > 0 for  q > c#,  k~'U2  is  decreasing  in  the 
direction  of  nereasing  f .ilong  the  line  a — a(P).  Since  k lU*  is  positive  and 
continuous  fc-  q > c*,  r~lU*  must  be  bounded  away  from  zero  for  q > ct. 
Hence  | k~l/tU  i is  bounc!«*d  away  from  zero  ona  ■ a{P),  0 > 0(P). 

In  terms  of  we  from  (3)  that 


(13) 


r'  V.  - *(kr'nu  + k'n*,) 


and  hence,  since  it  is  bounded  and  k — ♦ 0 for  a ■=  a(P),  <3  — » 0(P)  from  above, 
| k~u*it  | is  also  bounded  away  from  zero  for  a -»  a(P ),  0 > 0 (P),  0 sufficiently 
close  to  0 (P).  Therefore,  using  (4)  expanded  about  q = c *«  c„  we  have 


(14) 


Jk 


(?  ~ c.) 


T7i 


> 5 > 0 


for  a ■»  a(P),  0 > 0(P),  0 sufficiently  close  to  0(P). 

Now,  since  7 vanishes  at  P,  we  have,  by  (2),  i,(P)  ■»  0.  By  Taylor’s 
theorem,  since  the  second  derivatives  of  i exist  and  are  bounded, 


(15)  i,  « (9  - c.)F,  + (0  - <KP))F, 

where  P,  and  P,  are  bounded  functions  of  q and  0.  Along  the  characteristic 

a — a(P),  we  have,  by  (6) 


(16) 


0 - 0(P) 


f V?/C*  - 


1 q~l  dq  = m(q)(q  - c,) 


where  m is  bounded. 
(17) 

holds  and  therefore 


Hence,  for  a — o(P),  0 > 0(P) 
i . - (Pi  *f  mP,)(g  - cj 

— ° — is  bounded.  But  then 


9 “ e* 




(9  - e.),/4 


0 for 


a ■«  a(P),  0 — » 0(P)  and  this  contradicts  (14).  Thus  7 > 0 for  ? — c„  as  well 
as  for  q > c,  and  the  lemma  is  proved. 

It  remains  only  tc  apply  the  above  lemma  to  the  two  situations:  flow  past 
a fixed  airfoil  and  the  construction  of  airfoils.  Wc  treat  both  these  situations 
at  once,  and  note  first  that  the  boundary  of  the  domain  of  the  flow  in  the  hodo- 
graph  plane  is  the  streamline  corresponding  to  the  profile;  since  the  sonic  line 
intersects  the  profile  twice,  we  need  only  show  that  7 > 0 on  the  profile  to  prove 
that  7 > 0 for  <7  > c*. 

The  curvature  K of  any  streamline  is  given  by  the  formula 


K «=  d0/d»  =*  0.  cos  0 + 6,  8 in  0. 

But  the  transformation  x = x(q,0),  y = y(q,0 ) gives,  if  7 does  not  vanish, 

70.  = y . 

70,  = — xt 
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and  hence 

KJq  — y,q  cos  6 — x,q  sin  6 

- — - (y,+,  + x,t.)  “ - ~ > 

p p 

Lsing  (2)  we  then  find 

(18)  A'VV 

on  the  profile,  for  q > c,  > c. 

Consider  a set,  Af,  of  Mach  numbers  at  infinity  and  a set  of  profiles,  which 
depend  continuously  on  Af,  and  nave  finite  curvature.  Let  fC.  be  the  upper 
bound  for  the  curvatures,  and  suppose  that  J > 0 on  each  profile  for  Af  < Af  *; 
J — 0 for  Af  — Af*  at  seme  point  on  the  supersonic  section  of  the  correspond- 
ing profile.  Then,  by  (18) 

J > 1 /Klq  for  3/  < Af*, 

where  q is  the  limit  speed.  Therefore,  since  J is  a continuous  function  of  Af, 

J > 1 /Klq  >0  for  Af  - Af*; 

i.e.,  J cannot  vanish  on  the  boundary  streamline.  This  fact  is  derived  by 
Friedrichs  in  (3J. 

The  lemma  may  now  be  applied  to  show  that  J > 0 throughout  the  super- 
sonic region  and  on  the  sonic  line.  We  conclude  therefore  that  a limiting  line 
cannot  appear  in  the  flow  about  a fixed  airfoil  with  a bounded  supersonic  region 
adjacent  to  it,  and  if  it  appears  in  the  construction  of  profiles  with  adjacent  super- 
sonic regions  then  the  curvature  of  the  corresponding  profile  is  infinite  at  some 
point. 
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